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ABSTRACT 

In  this  paper,  the  existence  of  limit  cycles  for  a  class 
of  nonlinear  systems  is  explored.  Based  on  the  time- 
domain  approach  with  differential  and  integral 
inequalities,  the  phenomenon  of  the  stable  limit  cycle 
can  be  accurately  verified  for  such  nonlinear  systems. 
Furthermore,  the  exponentially  stable  limit  cycles, 
frequency  of  oscillation,  and  guaranteed  convergence 
rate  can  be  correctly  calculated.  Finally,  some 
numerical  simulations  are  provided  to  demonstrate  the 
feasibility  and  effectiveness  of  the  main  results. 

Keywords:  Limit  cycle,  nonlinear  systems,  stable  limit 
cycles,  exponential  convergence  rate 

1.  INTRODUCTION 

Nonlinear  network  may  cause  oscillations  with  fixed 
period  and  fixed  amplitude.  These  oscillations  are 
named  limit  cycles,  e.g.,  RLC  electrical  circuit  with  a 
nonlinear  resistor  and  Van  der  Pol  equation.  Limit 
cycles  are  special  phenomenon  of  nonlinear  networks 
and  have  been  widely  investigated;  see,  for  example, 
[1-12]  and  the  references  therein. 

Prediction  of  limit  cycles  is  very  meaningful,  in  view 
of  the  fact  that  limit  cycles  can  occur  in  any  kind  of 
physical  system.  Frequently,  a  limit  cycle  can  be 
worthwhile.  This  is  the  case  of  limit  cycles  in  the 
electronic  oscillators  utilized  in  factories  and 
laboratories.  There  are  at  least  four  methods  to 
explore  the  phenomenon  of  limit  cycles,  namely 
describing  function  technique,  Poincare-Bendixson 
theorem,  Piecewise-linearized  methodology,  and 
Lyapunov-like  approach.  The  disadvantages  of  the 
describing  function  method  are  related  to  its 


approximate  nature,  and  include  the  possibility  of 
inaccurate  predictions.  Besides,  the  Poincare- 
Bendixson  theorem  only  provides  a  necessary 
condition  to  ensure  the  existence  of  limit  cycles. 
Therefore,  even  the  conditions  of  the  Poincare- 
Bendixson  theorem  are  meted  for  a  system,  the 
existence  of  limit  cycles  cannot  be  guaranteed  for 
such  a  system. 

In  this  paper,  based  on  the  time-domain  approach  with 
differential  and  integral  inequalities,  the  phenomenon 
of  the  stable  limit  cycle  will  be  accurately  verified  for 
a  class  of  nonlinear  systems.  Meanwhile,  the 
exponentially  stable  limit  cycles,  frequency  of 
oscillation,  and  guaranteed  convergence  rate  will  be 
calculated.  At  last,  several  numerical  simulations  will 
offered  to  show  the  feasibility  and  effectiveness  of  the 
obtained  results. 

2.  PROBLEM  FORMULATION  AND  MAIN 
RESULTS 

In  this  paper,  we  consider  the  following  nonlinear 
system: 

X\(t)  =  -bx 3  (t)  -  X,  (?)[x  2  (t)  +  l] 


•  [xf  (; t )  +  x2  (t)  -  a]  , 

(la) 

II 

>< 

+ 

(lb) 

x3(t)  =  bxx  (t)  -  x3(?)[x22(7)  + 1] 

•  [x2  (t)  +  x2 (/)-«] , 

(lc) 
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x(t):=  [x,(t)  x2(t)  x3 (/ )]r  e'J? 3 


(Id) 


is  the  state  vector, 


x(o)r  :=  [x 

Where 

[X|°  XlJ>  is  the  initial  value,  and  a,b<E  'IJ 

represent  the  parameters  of  the  system,  with  a>  °. 
Clearly,  x  =  0  is  an  equivalent  point  of  system  (1),  i.e., 
the  solution  of  system  (1)  is  given  by  x^  =  0  if 
x(o) = 0  -pQ  avojq  {hg  apparent  case  of  x(°)  =  0  ,  in  the 
following,  we  only  investigate  the  system  (1)  in  case 
0f  *(o)  *  0 

Definition  1 


5(x)  =x\  +  x\ 


! .  Then  the  time  derivatives  of 


s  (x) 


and 


along  the  trajectories  of  system  (1)  is  given  by 

ds2(x{t )) 


dt 


-  =  2s(x(7))-  (2x1x1  +  2x3x3) 


=  -4  (xf  +  X3)(x2  +  l)s,2(x(/)). 

d6{x(t ))  _  x3Xj  -XjX3 


(2) 


dt 


2  2 
Xj  +  x3 


-  =  b. 


It  follows  that 

&(x(t))  -bt  +  tan-1 


r  \ 


(3) 


Consider  the  system  (1).  The  closed  and  bounded  In  the  following,  there  are  three  cases  to  discuss  the 
manifold  5(x)  =  0;  jn  the  x>~x3  plane,  is  said  to  be  an  trajectories  of  the  system  of(l). 
exponentially  stable  limit  cycle  if  there  exist  two 

positive  numbers  01  and  P  such  that  the  manifold  of  Case  1:  A  (°)  +  *3  (°)  =  «  (or  equivalently;  iS'M0)) = 0 ) 
5(x)  =  0  ajong  trajectories  of  system  (1)  meets  the 
following  inequality 


In  this  case,  from  (2),  it  can  be  obtained  that 

ds2(x(t)) 


|s(x(0)|  <  P  ■  exp[—  a(t-t0 )],  \/t>t0. 

In  this  case,  the  positive  number  a  is  called  the 
guaranteed  convergence  rate. 

Now,  we  are  in  a  position  to  present  the  main  results 
for  the  existence  of  limit  cycles  of  system  (1). 

Theorem  1. 

All  of  phase  trajectories  of  the  system  (1)  tend  to  the 


dt 


=  0, 


which  implies 


xf  (7)  +  x3  (7)  =  a,  V  t  >  0. 

Hence  we  conclude  that 

xft)  =  y[a  cos 


(4) 


x2(t)  =  -J~a  sin 


- 

/  A 

bt  +  tan-1 

X30 

\x\o)_ 

- 

f 

bt  +  tan-1 

X30 

lxio  )_ 

V  £  >  0, 


V  t>0. 


exponentially  stable  limit  cycle 


s(x)  =x2  +  x3  -  a  =  0 


m 


s(x(7))  =  0,  V  t  >  0, 

in  view  of  (3)  and  (4). 


the  1  3  plane,  with  the  guaranteed  convergence  rate 


,v(x(0))>0 


a:-< 


00, 

2a , 


2(^io  X30  )  tf 


tf  *10  X30  — 

if  x*  +  x30  >  a , 
xf0  +  x30  <  a. 


Besides,  the  states  X|  and  X]<^  exponentially  track, 
respectively,  the  trajectories 


y[a  cos 

bt  +  tan  1 

l* 

1° 

yfa  sin 

bt  +  tan  1 

L* 

1° 

^  X10  ) 

■  and 

\  xio  )_ 

,  in  the 

time  domain,  with  the  guaranteed  convergence  rate 

a/2 

Proof.  Define  a  smooth  manifold  5W  =  0  and  a 


continuous 


function 


6{x)  :=  tan  1  — 
x. 


with 


Case  2:  (0)  +  x* (0)  >  a  (or  equivalently;  s^» >  u ) 

In  this  case,  from  (2),  it  can  be  obtained  that  5  is 
a  strictly  decreasing  function  of  t  with 

/(x(0)>0,  V/>0  ^  and 

<  -4[x2(t)+x32(/)]s2(x(0) 

=  -4a  s2  (x(t)),  V  /  >  0. 

Applying  the  Bellman-Gronwall  inequality  with 
above  differential  inequality,  one  has 

s2(x(/))<s2(x(0))  •  exp[-  4 at],  V  t  >  0, 

This  implies 

K*(o)|  ^  Kx(o))|  •  exp[-  2  at\  V  /  >  0, 
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^xl{t)  +  x]{t)  -4a 

<  ^ X 2  0)  +  X3  (Y)  —  ^ X 2  (t')  +  X3  (t) 

=  |xf0)  +  X3  (0_^| 

=  |s(x(0)| 

<  |s(x(0))|  •  exp[-  2at\  V/>  0. 

It  yields 

x2  (0  +  x2  (t)  -  ra 

<  ^|s(x(0))|  -exp[-^],  Vt>0. 


+  u 


(5) 


Consequently,  by  (3)  and  (5),  we  conclude  that 

/ 

bt  +  tan-1 


JCj  (t)  -  yfa  cos 


{  \ 
X 


cos 


/  \ 


Z?/  +  tan  1 


Vxio  J 


—  yfa 


a  cos 


- 

( 

bt  +  tan-1 

X30 

l  X10  J  ] 

^(O  +  Xjh)  --Ja 

-Jxfo  +  Xjh)  -4a 
fiml  ■  exp[-  at\  V  t  >  0, 


r  w 

cos 

Z?/  +  tan  1 

X30 

^  X10  J  \ 

—  y[a 


a  sin 


(  w 

sin 

Z^  +  tan  1 

X20 

l  xio  J  J 

—  =  ~4[X<  (0  +  X3  (Of  2  (0+  if  (x(0) 

<  -4[x,2  (t)+ xj  (t)f  (x(t)) 

<  -4(x,20  +  x20  )s'2  (x(/)),  V  t  >  0. 


Applying  the  Bellman-Gronwall  inequality  with 
above  differential  inequality,  one  has 

s2(x(t)) 


<  s 


:  (x(0))  •  exp[-  4(x20  +  X20 }] ,  V  t  >  0, 


this  implies 

K*P))| 

<  |s(x(0))j  ■  exp[-  2(x20  +  x320 )] ,  V  t  >  0, 

2 

^ xl  (t)  +  x3  (?)  ~ 

^  X\  (0  4-  x 3  (0  —  yfci  ■  ^ x\  0)  +  X3  0)  +  yfu 

=  lxfOi  +  ^3  (0_a| 

=  Kx(o)| 

<  |s(x(0))|  •  exp[-  2(^20  +  x30  )}  V  t  >  0. 

It  yields 

ij xf  (t)  +  X l  ( t )  -  y[a 

VK*(°))|  •  exp[_  (*10  +  x30  )  J  V  t  >  0. 


(6) 


Consequently,  by  (3)  and  (6),  we  conclude  that 


yfa  sin 

bt  +  tan  1 

f  \ 

X20 

xx(t)- -Ja  cos 

Ztf +  tan  1 

/  W 

X30 

l  X10  )_ 

^X10  )\ 

/  \' 

(  \ 

tJ xf  (t)  +  x320)  sin 

Z^  +  tan  1 

X20 

l  X10  )_ 

— 

^X12(/)  +  X320)  cos 

a 

+ 

X30 

LX10  )_ 

bt  +  tan  1 

(  \ 
X20 

-  y[a  cos 

Z?/  +  tan  1 

/  A 

X30 

lxio  J _ 

lxio  J 

■j x2  (t )  +  x2  (?)  -  -J~a 

x2  (t)  +  x2  (?)  -  Va 
<  ^  •exp[-a/],  V^>0. 

Case  3:  4(0)  +  @(0)<«  (or  equivalently;  5W°))<0) 


(  ^1 

COS 

Z^  +  tan  1 

X30 

Lxio  J \ 

■Jxf(t)  +  x2  (?)  -  4a 
yjxt(t)  +  x*(t)-4a  | 

<  J|s(x(0))|  •  exp[-  (x2q  +  x20  >1  V  ?  >0, 


In  this  case,  from  (2),  it  can  be  obtained  that 


2(x(0) 


is 


a  strictly  decreasing  function  of  t  with 

*•(40)20,*,  >0>and 
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sin 


bt  +  tan 


/  \ 
X30 

n/3  COS 

^  n 

5 1 - 

-v/3  sin 

^  71 

5 1 - 

V  xio 

4 

and 

4 

,  in  the  time  domain, 


^ x i  (Y)  +  (/)  si 

-yfa  sin 


bt  +  tan  1 


V  10  y 


- 

(  ^1 

bt  +  tan-1 

X30 

l  X10  )  \ 

x j2  (^)  +  (^)  — 

xj2  (7)  +  X3  (/)  -  ^[a 


(  ^1 

sin 

Z^  +  tan  1 

X30 

l  xio  ) \ 

<  ^(xCO)!  •  exp[-  (4  +  Xjq  )}  v  t  >  0. 

This  completes  the  proof.  □ 


Remark  1.  It  should  be  pointed  out  that,  by  Theorem 
1,  the  system  (1)  can  be  regarded  as  nonlinear 

oscillators  with  the  amplitude  ^  and  the  frequency 
b.  Such  an  oscillation  is  entirely  independent  of  the 
initial  condition  and  limit  cycles  of  such  an  oscillation 
are  not  affected  by  parameter  variation. 


with  the  guaranteed  convergence  rate  2  .  Some 

state  trajectories  of  such  a  system  are  depicted  in 
Figure  3  and  Figure  4. 

4.  CONCLUSION 

In  this  paper,  the  existence  of  limit  cycles  for  a  class 
of  nonlinear  systems  has  been  considered.  Based  on 
the  time-domain  approach  with  differential  and 
integral  inequalities,  the  phenomenon  of  the  stable 
limit  cycle  can  be  accurately  verified  for  such 
nonlinear  systems.  The  exponentially  stable  limit 
cycles,  frequency  of  oscillation,  and  guaranteed 
convergence  rate  can  also  be  correctly  calculated. 
Finally,  some  numerical  simulations  have  been  given 
to  demonstrate  the  feasibility  and  effectiveness  of  the 
main  results. 
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3.  NUMERICAL  SIMULATIONS 


Example  1:  Consider  the  system  (1)  with  ^a,b^  ■2’ 4) 

and  x(°)  =  [2’0’2]  .  By  Theorem  1,  we  conclude  that 
the  phase  trajectories  of  such  a  system  tend  to  the 

exponentially  stable  limit  cycle  5W=xi  +  y  -2  =  0  -n 
the  F  plane,  with  the  guaranteed  convergence  rate 
a  =  4.  Besides,  the  states  F(0  and  F(0  exponentially 


track,  respectively,  the  trajectories 


V? 


cos 


4t  +  — 
4 


and 


■J2  sin 


in  the  time  domain,  with  the  guaranteed 


convergence  rate  2  .  Some  state  trajectories  of 

such  a  system  are  depicted  in  Figure  1  and  Figure  2. 


Example  2:  Consider  the  system  (1)  with  ^a,b- 

and  x(0)  =  [0.1,0, -O.l]  gy  Theorem  i;  we  conclude 
that  the  phase  trajectories  of  such  a  system  tend  to  the 

exponentially  stable  limit  cycle  ‘v(a4  =v,  +  x3  -  3  =  0  -n 
the  F  _F  plane,  with  the  guaranteed  convergence  rate 
« =  0.28  jn  addition,  the  states  F(0  and  F(0 
exponentially  track,  respectively,  the  trajectories 
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Figure  1:  Typical  state  trajectories  of  the  system  (1) 
withMM2,4)and*(0)  =  [2,0,2j\ 


Figure  2:  Typical  phase  trajectories  of  the  system  (1) 
with  M=  (2, 4)  and  A-C0)  =  [2, 0, 2f 


Figure  3:  Typical  state  trajectories  of  the  system  (1) 
with  M  =  (3,5)  and  x(0)  =  [0.1, 0,-0.  if 


Figure  4:  Typical  phase  trajectories  of  the  system  (1) 
with  M)  =  M  and  *(0)  =  [o.l,  0,-0.  if 
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